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Abstract. A new type of double-resonance spectroscopy of a quantum gas based on interaction-induced 
frequency modulation of a probe transition has been considered. Interstate interaction of multilevel atoms 
causes a coherence-dependent collisional shift of the transition between the atomic states |1) and |2) due 
to a nonzero population of the state |3). Thus, the frequency of the probe transition |1) — |2) experiences 
oscillations associated with the Rabi oscillations between the states |1) and |3) under continuous excitation 
of the drive resonance |1) — |3). Such a dynamic frequency shift leads to a change in the electromagnetic 
absorption at the probe frequency and, consequently, greatly enhances the sensitivity of double-resonance 
spectroscopy as compared to traditional "hole burning" , which is solely due to a decrease in the population 
of the initial state |1). In particular, it has been shown that the resonance linewidth is determined by the 
magnitude of the contact shift and the amplitude of the drive field and does not depend on the static field 
gradient. The calculated line shape and width agree with the low-temperature electron-nuclear double- 
resonance spectra of two-dimensional atomic hydrogen. 

PACS. 3 3.40.+f, 32.70.Jz, 34.50.Cx, 67.63.Gh, 



I 1 Introduction 

. Interaction of atoms and, more specifically, short-range 

■ or contact interaction, leads to a shift of internal atomic 
transitions commonly referred to as a contact shift. In the 

' case of multi-level atoms, the contact shift also depends 

■ on interstate coherence and populations of the states not 
. involved in the particular resonance. More specifically, ac- 
' cording to experiments on ^Li [1] and [2] and recent 
, theory [3l|4] , an incoherent population of such a third state 

■ causes a nonzero shift for fermions. On the other hand, 
I fermions in identical internal states do not interact via 
. s waves; therefore, a fully coherent ultra-cold Fermi gas 
' exhibits a contact shift only in the case of spatial inhomo- 
' geneity [SUS] . On the contrary, the respective contribution 

to the contact shift of bosons is the highest if the internal 
states are coherently populated and decreases by a factor 
of 2 as the sample decoheres. Thus, in a Bose gas, excita- 
tion of a drive transition to the state |3) induces frequency 
modulation of the probe transition 1 1) — 12) associated with 
Rabi oscillations between the states |1) and |3) [4]. The lat- 
ter leads to interaction-enhanced double resonance when 
two transitions are excited simultaneously. In this work, 
we consider this effect in the relaxation-free approxima- 
tion, analyze a possible double-resonance lineshape and 
compare it with experiments in spin-polarized atomic hy- 
drogen [7j£ . 

Interaction of multilevel systems with several resonance 
fields is commonly considered with the use of the evo- 



lution equations for the components of spin density ma- 
trix [9l. Inclusion of the population-dependent frequency 
shift makes these equations essentially nonlinear and their 
analytical solution impossible. Our approach is less com- 
prehensive but provides quantitative and physically clear 
description in the case of interest. 

2 Interaction-enhanced double resonance 

Physics of this novel phenomenon becomes transparent al- 
ready in the case of just three levels (Fig. [T]). In general, 
excitation of the probe transition |1) — |2) affects the fre- 
quency of the drive 1 1) — |3) resonance and vise versa. This 
may lead to specific effects and produce peculiar double- 
resonance spectra, which will be considered elsewhere. 

Here, for simplicity and clarity, we restrict ourselves 
to the case (see below for the exact condition) when the 
contact shift vanishes in the absence of the third state. 
Another advantage of this case is that it allows direct 
comparison with experiments on electron-nuclear double 
resonance (ENDOR) in atomic hydrogen [7j|8]. To simplify 
our consideration even further, we assume that the popu- 
lation of the state |2), which is the final state of the probe 
transition, is kept negligibly small and therefore the fre- 
quency of the drive transition remains constant. On the 
other hand, the frequency shift of the probe transition 
1) — 1 2) in the presence of the state |3) is 4 

hAcjf^li) = 2ni(A+ - An) + 2^2 (A22 - A+ ) 
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Fig. 1. (Left) Three- level scheme and (right) effect of Rabi 
oscillations between the states |1) and |3) on the |1) — |2) tran- 
sition frequency and intensity. Here, for better visibility, the 
Rabi period 2ii/fi\i is thought to be much longer than the 
time needed to detect the |1) — |2) resonance line. The oppo- 
site case is considered in Sees. |3] and |4] below. 

+ 2n3|C+3p(A+3-A+3), (1) 
HAojI^^^ = 2n3|Cr3|2(A23 - Ar3). (2) 

Here, rii is the population of the state A,^ = {ij\\\ij)± 
is the matrix element (excluding the spatial factor) of 
the interaction strength A — ATrh?a/m commonly used 
in a cold-collision regime, m is the atomic mass, and a 
is the respective scattering length. Cf^ are the normal- 
ized (|Cj^P -I- \C~j\'^ = 1) amplitudes of the {+) symmet- 
ric and (— ) antisymmetric component of properly sym- 
metrized diatomic wavefunction |ki,qj) = C^ip±\ij)+ + 

C~j'4>zf\ij)-, where the upper (lower) subscript stands for 
bosons (fermions) and the internal (pseudo-spin) and spa- 
tial parts are, respectively, \ij)± = "Tj'-I*-'^ l-'*^-' ^^'^ 
V'± = ;75(V'k(ri)V'q(r2) ± V'q(ri)V'k(r2)). 

Equation ([2]) should be compared with Eq.(l) of Gupta 
et al. [T] and Eq. (1) of Regal and Jin 2 , which hold in 
a fully incoherent case of |C{3p = ^, as well as with the 
general theoretical result (Eqs. (8)-(9) of Baym et al. [3]). 
Equation ([T]) generalizes the two-level formula (4) of Gib- 
ble [5]. Clearly, the shift in the two- level gas (713 = 0) is 
zero for fermions and for bosons with An = A^2 = •^22- 
The latter is exactly the case mentioned above and agrees 
with our previous result [TU] and with the observations 
in atomic hydrogen [71lll|. Moreover, in agreement with 
other authors pilSl lT^IT^ . the shift is independent of the 
sample coherence, because the fraction |C]^| of pseudo- 
singlets does not enter Eqs. (P) and Thus, we come 
to the well-known conclusion that the interstate coher- 
ence in a spatially homogeneous two-level gas cannot be 
probed by the contact shift [5]. Note, that spatial inho- 
mogeneity renders the atoms distinguishable and causes 
specific coherence-dependent shifts of opposite signs for 
bosons and fermions [SJIB]. 

The situation in a three-level gas (713 7^ 0) is funda- 
mentally different. We see that the interstate coherence of 
a spatially homogeneous gas can be probed by the contact 



shift in the presence the third state, like it has been done 
in fcrmionic ^Li [T. It is also remarkable that a homoge- 
neous Fermi gas does exhibit a nonzero clock shift if the 
two states involved in the transition are not fully coher- 
ent with the third one (C{^ 7^ 0), the latter is populated 
(713 ^ 0) and the respective scattering lengths are differ- 
ent (aj^3 7^ 023)- this context, an incoherently populated 
third state acts as a buffer gas of foreign atoms. One can 
also see that, like in the case of just two levels, the co- 
herence between the one-atom states |1) and |2) coupled 
by the resonance transition does not affect the clock shift 
Z\aji2(3) of this resonance. 

In the absence of relaxation, the populations of the 
states |1) and |3) evolve (see Fig. [T]) under continuous ex- 
citation of the transition |1) — |3) in the initially |l)-state 

sample as = nsin^ ^13 sin^ ni — n ^ n^, where n 
is the total density of the gas, sinf?i3 = fli^/fiis, J?i3 = 
y/ Aujf^ + ^2^3, f2i3 is the respective Rabi frequency, and 
Z\cji3 = Wd — 1^13 is the detuning between the resonance 
frequency W13 and the frequency ojd of the drive field. 
Note, that the gas remains fully coherent, since pseudo- 
spins 5*13 of all atoms experience coherent precession on 
the |1) — |3) Bloch sphere around the tilting angle 6*13 and 
stay parallel to each other (Fig. [1]). Thus, the resonance 
frequency CJ13 remains unchanged while the |1) — |2) tran- 
sition is dynamically shifted to the frequency 0] 



huji2 — huj'^ + 2nZiAsin'^ 6*13 sin^ — (3) 



where AX = AJ3 — X^^. The amplitude of this frequency 
modulation can easily be comparable with or even much 
greater than the linewidth of the probe transition. Thus, 
excitation of the transition |1) — [3) periodically drives 
the probe transition out of resonance. For clarity, Fig. [T] 
corresponds to slow driving in a sense that the Rabi pe- 
riod 27r/f2i3 of the drive transition is thought to be much 
longer than the time needed to detect the |1) — |2) reso- 
nance line. The opposite case of fast driving is considered 
in Secs.[3]and|l]below. At slow driving, the entire |1) — |2) 
spectrum moves periodically fore and back along the fre- 
quency axis and simultaneously changes in amplitude. The 
waveforms shown in Fig. [1] are simply the snapshots of 
the spectrum at different phases of the Rabi cycle. Obvi- 
ously, electromagnetic absorption at the fixed frequency 
Wp of the probe field also changes periodically with a pe- 
riod of the |1) — |3) Rabi oscillations. It should be empha- 
sized that these changes are associated with changes in 
both the population of the initial state and the transition 
frequency, in contrast to conventional double resonance, 
which is solely due to a change in the population of the 
initial state caused by the drive transitions. Clearly, such 
interaction-induced frequency modulation can greatly en- 
hance the effect, which therefore can be called INteraction- 
Enhanced DOuble Resonance (INEDOR). 
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Field detuning (co -co )/y H 

Fig. 2. Field dependence of the [1) — |2) frequency shift (in 
units of ypHd) under the CW excitation of the |1) — |3) reso- 
nance. Horizontal axis is the field detuning (cjis — a;d)/7d from 
the |1) — |3) resonance in units of the excitation field H^. Solid 
and dashed line are, respectively, the upper bound of the sum 
of the Zeeman and mean-field contributions and the Zeeman 
contribution alone. Horizontal lines indicate (dash-dotted line) 
the value of cji2 at minimum and (solid line) the probe fre- 
quency. Vertical dashed line corresponds to the resonance field 
for the |1) — |3) transition. The parameters correspond to 2D 
atomic hydrogen with a density of 3 • 10^^ cm~^ in a high po- 
larizing field of 45 kG, except for the sign of the contact shift: 
Hd = 1 mG, AHc — 89 G, 7p = 7e, 7d = 7p- Inset is the 
overview of the same dependence. 

3 INEDOR spectrum 

Let us consider in more detail a possible line shape of 
the INEDOR spectrum. In contrast to the case of slow 
driving illustrated in Fig. 1, we assume that simultane- 
ous frequency and amplitude modulation of the |1) — |2) 
absorption line is relatively fast (which is the case in a suf- 
ficiently high drive field) and therefore integrated by the 
detection system. This implies that the |1) — |3) Rabi fre- 
quency is much higher than the rate of field or frequency 
sweep through the |1) — |2) resonance or the inverse time 
constant of the detection system, ilisT ^ 1. In prac- 
tice, as discussed in Sec. |4] below, a less severe condition 
^i3T 3> 1 is sufficient. Thus, we deal with time-average 
absorption at the probe frequency in a sense that the ab- 
sorption signal is integrated over many modulation (Rabi) 
cycles. 

Generally, the energies of all three levels and, conse- 
quently, both transition frequencies depend on the exter- 
nal static field. The nature of this field is insignificant but 
for definiteness we shall consider the magnetic field H. 
Since our analysis is conducted in terms of pseudo-spins, 
the respective contribution to the frequencies is the Zee- 



man splitting jd{p)H, where 7d(p) is the effective gyromag- 
netic ratio of the drive (probe) transition. Let the field Hq 
correspond to the exact |I) — |3) resonance, wi3(7?o) — <^d- 
Then the deviation h = H — Hq of the field from this value 
results in the drive frequency detuning AuJi3{h) = j^h. 
Owing to this Zeeman contribution to the effective pre- 
cession frequency J7i3, the amplitude of the oscillating 
component of wi2 in Eq. (jI2[) . i.e., the amplitude of the 
frequency modulation of the probe transition turns out be 
a Lorentzian function of the field 

^ f2nAX\ Hi 

where is the amplitude of the drive field Hd{t) = 
-ffd exp(iajdi), ^?i3 — 7dffd- On the other hand, the static 
(zero-density) component of uji2 also includes the Zeeman 
term 

i^^^ {Ho + h)= oj[l^ (Ho) + 7p/i. (5) 

The sum of these two contributions is an upper bound of 
the probe frequency, 

As a result, uji2(t) oscillates between the Zeeman-only 
lower bound (Eq. dS]), dashed line in Fig. ^) and the 
Zeeman plus mean- field upper bound (Eq. ^ , thick solid 
line in Fig. ^) at the field-dependent frequency Qij,{h) — 
jd\/Hj+h?. As is seen from Eq. the upper bound of 
the probe frequency is generally a nonmonotonic function 
of the static field. 

As mentioned above, the Rabi frequency J7i3 = i7i3(0) 
is assumed to be sufficiently high, Hi^t I, in which 
case the absorption amplitude is integrated by the de- 
tection system over many Rabi cycles. The time-average 
absorption amplitude A{h,LUp) within the bounds ^ and 
^ is proportional to the average population of the ini- 
tial state multiplied by the probability density p{h,iOp) — 
[Sli3{h)/2TT][duJi2{h,ujp)/dt]~^ to find the system at given 
h and Wp. It is readily shown that 

1 ~ X 

A{h,ujp) (X , (7) 

2t:\/ a;(sin^ disih) — x) 

where 

n[ujp - J^^ {Hp + h)\ 

and 1 — X are the relative populations of the states |3) and 
I), respectively. Thus, the probability density peaks at the 
lower {x — 0) and upper (x = sin^0i3(/i)) bound of the 
probe frequency. Note, that in the latter case absorption is 
further reduced due to a smaller population of the initial 
state and vanishes at x = 1 (which implies h — 0)^ i.e., at 
the maximum of the Lorentzian function (|4l). 

The spectrum that appears when both the drive and 
probe frequency are fixed and the field is slowly swept 
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through the resonance consists of up to four peaks corre- 
sponding to the field values, at which the probe frequency 
shown by the solid horizontal line in Fig. [5] crosses the 
upper (Zeeman plus mean-field) and lower (Zeeman only) 
bound of uji2 (points A, B, C, and D). The actual number 
of lines depends on the amplitude of the contact shift, the 
ratio 7d/7p of the Zeeman quanta of the excitation fields 
and the relation between the drive and probe frequency. 
The peaks lying within the drive resonance (i.e., within 
the segment AB in Fig. [2]) are much weaker because the 
atoms spend only a small fraction of time in the resonance 
conditions. This implies that the |1) — |2) absorption am- 
plitude is strongly decreased, which produces a "hole" in 
the resonance curve. 

So far we considered a spatially homogeneous system. 
There could be three types of inhomogeneity associated 
with the inhomogeneous (i) static field, (ii) excitation fields 
and (iii) gas density. In terms of the upper bound of the 
probe frequency shift in Fig. [21 these correspond to a 
spread in (i) the overall vertical position of the curve 
(equivalent to the inhomogeneity of the probe frequency), 
(ii) the width of the drive resonance and (iii) the ampli- 
tude of the drive resonance plus the overall vertical shift. 

To illustrate the effect of inhomogeneity we consider in 
more detail a linear gradient VH of the static field in an 
otherwise spatially homogeneous infinite sample. In this 
case, all possible field values are simultaneously present in 
different parts of the sample. Thus, the absorption ampli- 
tude at a given probe frequency Wp is proportional to the 
integral /(wp) = / A(/i, Wp) (|^) d/i of Eq. (O along the 
ujp = const line within the bounds of the probe frequency 
(i.e., over the horizontal segments AB and CD in Fig. 
Owing to a constant field gradient, integration over the 
field is equivalent to integration over space provided that 
the number of particles per unit volume n is replaced by 
the number of particles per unit field ^ oc n\\7H\^^. The 
result of numerical integration with the parameters of ex- 
periments with 2D atomic hydrogen [71I51IT5] is shown in 
Fig. [3] as a function of the drive frequency Wd (lower hor- 
izontal axis) at constant Wp. Alternatively, the INEDOR 
spectrum can be detected by sweeping Wp (upper horizon- 
tal axis) at constant Wd- In the latter case, the spectrum 
is inverted on the frequency scale because an increase in 
Wd corresponds to an increase in the resonance field Wd/7d 
and, consequently, to a positive displacement of the mean- 
field Lorentzian peak on the cjuih) curve (O in Fig. [5J 
The latter is equivalent to a negative change in the probe 
frequency ujp at constant Wd- Which way of observing the 
INEDOR spectrum is preferred depends on the details of 
a particular experiment. 

The absorption amplitude as a function of Wd (Fig. [3]) 
is nearly constant well below and above the drive reso- 
nance, decreases substantially within the resonance and 
has a sharp maximum at Wp = Wmin, i-e., when the mini- 
mum value of Eq. ^ coincides with the probe frequency. 
This explains the dispersion-looking ENDOR spectra of 
2D atomic hydrogen [7l|8] (see Sec|4] below). Physically, 
the hole in the absorption amplitude is because the atoms 
of the otherwise resonant region of the sample are dynam- 



ically driven out of the probe resonance and spend only 
a small fraction of time in the resonance conditions. As 
already mentioned, the probability of finding them in res- 
onance peaks at the lower ^ and upper (O bound of the 
probe frequency and, according to Eq. ([7]), is the lower 
the higher is the difference (HI between the two bounds of 
UJ12. In this respect, there is no surprise that the positions 
of the / minimum and the contact shift Q maximum 

nearly coincide, ujp ~ ujf2{HQ). On the other hand, the 
absorption maximum is due to the fact that the drive res- 
onance introduces zero gradient of the |1) — |2) transition 
frequency in a certain region of the sample. As a result, 
much more atoms become resonant. A decrease in the pop- 
ulation of the initial state |1) at the point of zero gradient 
turns out to be of minor importance for the chosen values 
of the parameters. However, according to Eq. ([6]), it can 
be significant at |2nZ\A| ~ h'ypHd, which occurs at low gas 
density n and strong drive field H^. 

The dispersion-looking INEDOR spectrum of a gas in 
a spatially inhomogeneous static field can be also qualita- 
tively explained in a different way. In the case of a positive 
contact shift, the parts of the inhomogeneously broadened 
probe resonance line are shifted towards higher fields by 
the amount proportional to the local amplitude of drive 
resonance curve. The shift increases with an increase in 
field at the low-field side of the drive resonance and de- 
creases at its high-field side. Thus, the low-field side of 
the spectrum is "rarefied" whereas the high-field side is 
"compressed" in a sense that the number of atoms per 
unit field, which are involved in the probe transitions at 
lower (higher) fields, becomes smaller (larger). Accord- 
ingly, the probe absorption intensity is decreased at the 
low-field side of the drive resonance and increased at its 
high-field side. When uj^ is swept from below at constant 
Wp, the high-field side of drive resonance comes first. That 
is, the probe absorption intensity is increased at lower and 
decreased at higher Wd, as illustrated by Fig. [31 

As discussed above, the width of the double-resonance 
curve in the probe-frequency units can be estimated as 
the difference cjmin —u!^{Hq) (see Fig. [2). For |2nZ\A| 
h'jpHfi the minimum of the right-hand side of Eq. (|6]) oc- 
curs at 

{Hj + _ 4nA\ 
Hjh hjp ■ 

Typically, h :$> and Eq. (© yields 

.3 . l!^. (9) 

%P 

This corresponds to 

duji2 = ujuih) -uj^"^{Ha) ~ --jph (10) 

and we finally arrive at 

Sio.s = ^Slui2 ^ lid{2AH,Hiy/^ (11) 
7p ^ 

where AHc — 2nZ\A(fi7p)^^ is the maximum contact shift 
of the |1) — |2) resonance in field units. Note, that the res- 
onance linewidth depends on the gas density and drive 
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Probe frequency [co - co *"'(// )]/27r, kHz 
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the \a) — \b) resonance. According to Eq. ([T|), the \b) — |c) 
resonance is dynamically shifted to the frequency |4J 
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Drive frequency {(a-ca^^yin, Hz 



800 



Fig. 3. |1) — |2) absorption amplitude as a function of drive 
frequency (expressed as the detuning, in Hertz, from the |1) — 
|3) resonance) for the delta-functional spectrum of the probe 
field. Sharp peak corresponds to the frequency minimum in 
Fig. [21 Parameters are the same as in Fig. [J] 



field as n^/^ij^/a j^gj^arkably, the width of the INEDOR 
spectrum is independent of the static field gradient. On 
the other hand, the spectrum intensity is inversely pro- 
portional to |ViJ|. A one-dimensional inhomogeneity of 
the drive field in the direction of \VH\ also does not 
affect the linewidth because in this case is a single- 
valued function of h. The effect of such an inhomogene- 
ity is that the mean-field peak ^ becomes asymmetric 
and the INEDOR linewidth is determined by the local 
value at the minimum of uji2 in Eq. ([5]). It is also worth 
mentioning that the linearity of the Zeeman contributions 
to aJi2 and wia is unnecessary because these terms can 
be always linearized in the vicinity of the resonance field. 
In this case, the gyromagnetic ratios should be replaced 
by the respective field derivatives, 7d dui^/dH and 
7p duu/dH. 



4 ENDOR in atomic hydrogen 

The spin states of a single hydrogen atom in a high mag- 
netic field are commonly referred to as \a) w I ^ +)i 

|&) = I ), |c) w I -f and \d) = \ + +), in the order of 

increasing energy (± signs stand for the projections of the 
electron and nuclear spin on the field direction; we also ne- 
glect small impurities of the opposite spin projections in 
the states \a) and |c)). Ahokas et al. [71[S] observed a nar- 
row feature in the ESR \b) — \c) line while sweeping through 



■ WbcO 



ml 2 



where n is the 2D density of adsorbed H atoms and I ~ 
0.5 nm is their out-of-plane delocalization length in the 
adsorption potential so that n/l stands in place of the 3D 
density in Eqs. (P; Aa = as - at = -30(10) pm [llfTl] is 
the experimental difference between the singlet and triplet 
scattering length of two ground-state hydrogen atoms (to 
be compared with the theoretical values ranging from —42 
to —55 pm [TSllTClfTT] ). The amplitude of the respective 
modulation of the resonance field in 2D atomic hydrogen 
at zero detuning {9ab = f ) is very large, AHc = 1.5(5) ■ 
10~^*^ G-cm^. At a typical 2D density of 3 • 10^^ cm'^, 
which corresponds to a 3D density of 6 • 10^^ cm~^, one 
has AHc ~ 90 G. The amplitude of the drive RF field is 
Hd ^ 10-3 G and we find from Eqs. ([T^ . and ([TT1) 
h ^ 5.7 • 10"^ G and Sujab/'^T^ ^ 350 Hz, very close to a 
numerical result of 330 Hz for the maximum-to-minimum 
distance in Fig. [31 This agrees fairly well (and, remarkably, 
without any fitting parameters) with the 120- Hz experi- 
mental width of the ENDOR absorption signal reported 
in [71IH]. One should take into account that the contact 
shift in atomic hydrogen is negative and therefore the ex- 
perimental spectrum is horizontally inverted with respect 
to the one shown in Fig. [31 

The experimental amplitude of the drive field corre- 
sponds to moderate Rabi frequency J7i3 ~ 30 . How- 
ever, the effective precession frequency i7i3 is two-three 
orders of magnitude higher in the most significant part 
of the spectrum (h/H^ ~ 560) except for a very narrow 
region near zero detuning. On the other hand, the fre- 
quency was typically swept through the ENDOR line for 
T ~ 1 s [T84. Thus, the condition J^isr » 1 of fast driving 
was fulfilled. Note, that the experimental spectrum can be 
affected by inhomogeneity of the drive field and gas den- 
sity across the 2D sample, as mentioned above. According 
to Eq. ([Hj), both lead to inhomogeneity in the position of 
the sharp peak in Fig. [31 which is therefore substantially 
broadened and reduced in amplitude. Actually, in experi- 
ments [S' , the amplitude of drive RF field was not exactly 
known and varied within the sample by a factor of three 
or so. 

It can be shown that relaxation and other possible in- 
elastic processes disregarded in the above consideration 
are of minor importance in atomic hydrogen [4J. Neither 
should the spectral width ^Wp/wp ^ 10^^ of the mm-wave 
source [5] produce any noticeable effect, as it is much nar- 
rower than the INEDOR linewidth 7p/i/wp ~ 3 • 10"^ in 
probe-frequency units. 

Ahokas et al. [71[H] also observed double resonance in a 
different way, by sweeping the static field at constant drive 
and probe frequencies and detecting the \h) — \c) transi- 
tion. This resulted in a change in the overall intensity 
of the absorption line depending on the drive frequency 
rather than in a narrow feature in the resonance curve, 
as one would expect from the above analysis. This issue 
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remains unclear to us. Here, we can only mention that the 
case i7^3 < 7d|9i?o/9i| of slow driving or, equivalently, of 
fast sweeping, which could take place in such experiments, 
should be analyzed separately. 

5 Conclusions 

Thus, we have shown that simultaneous excitation of two 
transitions in a cold gas of multi-level atoms can result 
in nonlinear double resonance induced by a change in the 
contact shift of the probe transition frequency due to Rabi 
oscillations between the states involved in the drive tran- 
sition. The linewidth of such interaction-enhanced double 
resonance is determined by the amplitude of the contact 
shift and the magnitude of the drive field and does not 
depend on the field gradient. The INEDOR line shape 
and width agree with those observed in 2D atomic hy- 
drogen. Ahokas et al. [71[5] speculated that the double- 
resonance spectra they observed in atomic hydrogen could 
be due to electromagnetically-induced transparency or co- 
herent population trapping, the well-known nonlinear ef- 
fects in three-level systems. The detailed discussion of 
this opportunity will be presented elsewhere. Here, we 
have shown that there is at least another realistic sce- 
nario based on interaction-induced double resonance. In 
a way, it can be also regarded as electromagnetically- 
induced transparency, although the physical mechanisms 
behind these two coherent phenomena are totally differ- 
ent. 

Reliable detection of conventional double resonance re- 
quires a significant (of order unity) depopulation of the 
initial state. As already explained in Sec. [I] INEDOR can 
be observed at a much lower depopulation because the 
probe transition frequency is driven completely out of res- 
onance by even a very small population of the state |3). 
More specifically, the amplitude of the frequency modula- 
tion should be higher than the homogeneous linewidth of 
the probe transition or the spectral width of the source, 
whichever is greater. For example, in the case of 2D atomic 
hydrogen, the latter is dujp/ujp ^ 10~^ [5], which corre- 
sponds to 2 • 10^ cm~^ or n^/n ~ 10~^. Owing to 
such an ultimate sensitivity, the specific double resonance 
technique proposed in this work could be an effective tool 
for studying interactions in atomic systems and beyond. 
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